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^ ' We analyze the generation of polarization cross-talk in Stokes polarimeters 

by atmospheric seeing, and its effects on the noise statistics of spectro- 
polarimetric measurements for both single-beam and dual-beam instruments. 
We investigate the time evolution of seeing-induced correlations between 
Q_i' different states of one modulation cycle, and compare the response to these 

correlations of two popular polarization modulation schemes in a dual-beam 
system. Extension of the formalism to encompass an arbitrary number of 
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modulation cycles enables us to compare our results with earlier work. 
Even though we discuss examples pertinent to solar physics, the general 
treatment of the subject and its fundamental results might be useful to a 



wider community. © 2011 Optical Society of America 
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1. Introduction 



Many important solar phenomena are driven by the interaction of turbulent plasma with 
magnetic fields. Measurements of the full state of polarization of solar spectral lines are 
routinely used to infer properties of solar vector magnetic fields from their imprints in the 
emergent spectra through the Zeeman and Hanle effects. All ground-based measurements 
are detrimentally affected by image distortions introduced by atmospheric seeing, which 
has significant power at frequencies up to at least 100 Hz, significantly higher than the 
frame rates of typical CCD cameras. The seeing-induced errors are expected to increase as 
telescope apertures increase in size: the telescope can resolve more structure (and hence 
steeper gradients) in images; the number of isoplanatic patches to be corrected through 
adaptive optical techniques increases with the square of the telescope aperture. 

The Advanced Technology Solar Telescope (ATST), a 4-m aperture, off-axis telescope, 
will soon enter construction phase. The European Solar Telescope (EST) is a 4-m class 
telescope currently under study. Spectropolarimetry is the prime mode of operation of this 



new generation of solar telescopes, but there is an urgent need to identify potential pitfalls 
before committing to a particular design for the modulation and detection of polarized light 
using such telescopes. 

Space-based spectro-polarimetric instruments, like those onboard the Hinode (Solar-B) 
and Solar Dynamics Observatory spacecraft, are by necessity fed by smaller telescopes, and 
while their observations profit greatly from the absence of atmospheric disturbances, they 
are still affected by residual image motion due to spaceship jitter. 

The purpose of the present paper is to re-examine and extend earlier theoretical studies 
of errors introduced by the effects of atmospheric seeing in order to help the design of polar- 
ization systems for these large, multi-national projects, and for future solar space missions. 

Two modes of operation are of particular interest. Slit-based spectro-polarimeters often use 
long integration times, including a number of modulation cycles that is typically of the order 
of a few tens. Tunable imaging spectro-polarimeters use instead much fewer measurements 
- often just one - of the modulated intensities, for each wavelength position. The formalism 
developed in this paper is general, and can be applied to both types of instruments and 
observations. 

The plan of this paper is the following. In the next section we lay the foundation of the 
formalism. Using the point of view of the statistics of random processes, in Sect. [2] we define 
the fundamental statistical observables for spectro-polarimetry. In Sect. El this formalism is 
extended to the case of dual-beam polarimetry, emphasizing those aspects of the problem 
where this extension is not trivial. In Sect. H] we make use of the conceptual framework 
of the statistics of stationary random processes to gain an insight on the effects of seeing 
correlations on spectro-polarimetric measurements. In particular, in that section we study 
qualitatively the effect of AO-type corrections on the temporal decay of seeing correlations. 
Those results are successively applied to study the performance of two popular modulation 
schemes in the presence of seeing. In Sect. [5l these results are further generalized to the case 
in which seeing correlations extend over the duration of one elemental observation (e.g., one 
slit position), and the effects of seeing on the performance of a balanced modulation scheme 
are studied both in the absence and in the presence of low-order (i.e., tip-tilt) AO correction 
(based on observed data). Finally, we compare our results and conclusions with those from 
previous works on the subject. 

2. Statistical description of the effects of atmospheric seeing 

We conventionally describe the polarization state of a radiation beam by its Stokes vector, 
S = (S±, S2, S3, S4), where S± is the intensity of radiation, S2 and S3 describe the two inde- 
pendent states of linear polarization, and finally S4 is the parameter for circular polarization. 
Because photon detectors are practically sensitive only to the intensity of the incoming 



radiation, the measurement of polarized radiation requires that its states of polarization be 
encoded into intensity signals. This encoding process is conventionally called modulation, 
and it is most often achieved through time-varying optical devices (modulators) that modify 
in a known way the polarization states of the incoming radiation. 

Let m(t) be the four-vector describing the polarization modulation process for the four 
Stokes parameters, so that m(t) ■ S gives the detected signal of the modulated intensity at 
time t (see Eq. [2] below). We will assume that this modulation vector is perfectly known 
(either by design or calibration), and that it is a periodic function of time, so that m(t + T) = 
m(t), where r is the period of the modulation cycle. Because of these properties, the time 
process of polarization modulation is both deterministic and stationary in a statistical sense. 

It is evident that one needs at least four independent measurements during a modulation 
cycle to fully determine the Stokes vector. The actual number of measurements n taken 
during a given modulation cycle represents the number of modulation states for that cycle. 

Let S be the incoming Stokes vector averaged over the spatial resolution element in the 
absence of seeing distortion. We assume that the physical conditions of the light emitting 
region do not vary over the time interval T needed to perform a spectro-polarimetric ob- 
servation of a given element of spatial or spectral sampling with the required sensitivity^] 
Therefore S is a function of the pixel position over the observed area that is assumed not to 
vary during this characteristic time. 

In the presence of atmospheric seeing, a given pixel on the detector collects photons from 
different elements of the observed region, with a characteristic time scale to (~ 0.01 s). The 
spatial distribution of the Stokes signal coming from these different elements is, to lowest 
order, completely determined by the true Stokes vector S that would be detected by a given 
pixel in the absence of seeing, and by its gradient VS. Therefore, we can assume that, at 
any given time within the duration T of a particular observation, the Stokes vector on the 
observed region that is mapped to a given pixel on the detector is expressed by [H[2] 

S(t) = S + VS • x(t) , (1) 

where x(t) is the average displacement vector on the plane of the skjo (e.g., expressed in 
fraction of the pixel unit) due to atmospheric seeing, whereas by assumption S and VS are 
constant during the observation. 

If only one isoplanatic patch contributes to the detected signal at any given time, then the 



1 In what follows, the term "observation" will always refer to such elemental operation, for instance, the 
observation of one slit position of a spatial scan, in the case of a grating-based spectro-polarimeter, or of one 
wavelength position of a spectral scan, in the case of a tunable, imaging spectro-polarimeter. 

2 In this study we ignore defocus by seeing, for which a full treatment of the diffractive properties of the 
entire optical system atmosphere+telescope would be necessary. Thus, here we are only concerned about 
lateral displacements of the image point by atmospheric turbulence. 



displacement vector x(t) can be efficiently compensated through a low-order, tip-tilt correc- 
tion. In the case of large telescopes, instead, a numbeio of different isoplanatic patches will 
contribute to the detected signal, so that the displacement vector x(t) in Eq. ([T]) must be 
regarded as an average of the various displacement vectors associated with the different con- 
tributing patches. In this case, an effective compensation of the seeing can only be achieved 
by use of high-order adaptive optics. However, S and VS remain the same for all patches, 
and therefore the model of Eq. (TjQ) still applies in this more general case. 

The modulated intensity recorded by an ideal detector (i.e., neglecting bias and read-out 
noise), for the i-th state of the modulation cycle, centered around the time £;, and integrated 
over the exposure time At, is given by 

4 +At/2 

J i = K ^2 I mj(t + t t ) [Sj + VSj ■ x(t + U)] dt + 5Ji, i = 1, . . . , n , (2) 

j=1 -At/2 

where k is a dimensional scaling constant for the detector, n is the number of modulation 
states in the cycle, and 8^ is a random fluctuation due to photon noise statistics. 

It is reasonable to assume that seeing can be regarded as a stationary and ergodiqj random 
process [I]. It therefore makes sense to derive the expressions for the expectation value and 
variance of ^ - respectively, E(j?i) and a 2 (J?i) - as these provide important information 
about the quality of the observations. 

Because of the random nature of atmospheric seeing, E(x(t)) = 0. Similarly, the photon 
noise fluctuation of the measurement has zero expectation value. Therefore, using Eq. (J2J), 
and the fact that the modulation process is assumed to be fully deterministic and stationary, 
we obtain 

a r + Ai / 2 

mtjSj + ±_ J mi{t + u) vs r EW + ti ))* 

L -At/2 



A = Ei^) = K AtJ2 



nAt y^ rriij Sj , i — 1, . . . , n , (3) 



where we have defined 



+A4/2 

miJ ~At / m i( t + t i) dt - ( 4 ) 

-At/2 



3 This number is of the order of (D/ro) 2 , where D is the telescope's diameter, and r® is the Fried parameter 
of the atmospheric seeing [3] . 

4 Simply stated, ergodicity is the statistical property of a dynamical system that allows to use time averages 
in place of ensemble averages. In other words, it is characteristic of an ergodic process to attain, during its 
temporal evolution, all of its possible configurations, with a probability distribution identical to that of the 
ensemble average. 
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The nx4 matrix M = (m^) so constructed is conventionally called the modulation matrix. 

Equation ([3]) shows that the expectation value of the i-th modulated intensity signal is 
determined exclusively by the physical Stokes vector S, regardless of the modulation scheme 
adopted. It is important to note that, in both Eqs. (J3J) and (J3J, the subscript i no longer 
represents a specific instant ti, like in Eq. fl2]), but strictly only the corresponding step position 
of the n-state modulation cycle. In particular, for Eq. (jlj), this is a direct consequence of the 
stationarity of the modulation process. 

Similarly, the variance of ^, for % = 1, . . . , n, is given by 



4 +A '/ 2 1 


2 


k^2 / rrij{t + ti) VSj -x(t + ti) dt + SJ^i 




j=1 -At/2 





4 +At/2 

= k 2 Y^ // m i(* + *i) "**(*' + **) ^([ V5 i • x (t + ^)][ v ^fc • «(*' + ti)}) dtdt' 
j,k=i 



At/2 



+<#•*) , (5) 



where we indicated with a p (j^) the RMS photon noise. 

Equation (jSJ) has no immediate applicability to Stokes data analysis, since it relies on 
physical parameters, such as V5, and x(t), which are not directly available from typical ob- 
servations. Instead, the quantities on the l.h.s. of Eq. (|5]) are determined in practice from the 
statistics of repeated intensity measurements over many modulation cycles p. Nonetheless, 
Eq. ((SJ) helps clarifying the physical origin of the noise in the measurement of the modulated 
intensities, and therefore can provide insights on optimal choices of modulation schemes and 
frequencies for reducing the final error. 

We note that also in Eq. ([5]) the index i no longer refers to a specific instant in time, ti. 
Once again, this follows from the periodicity of m(t), and the stationarity of the random 
process described by x(t). In fact, it is possible to drop altogether the temporal shift ti in the 
arguments of x(t) above, since the expectation value in the second line of Eq. (jSJ) depends 
only on the difference of those two arguments, if x(t) is stationary (see Sect. HI). 

Following [5], the optimal demodulation matrix (under specific assumptions; see comments 
at the end of this section), which allows us to infer the incoming Stokes vector from the 
modulated intensity signals, is given by D = (dij) = (M'M) _1 M*. If we indicate the measured 
Stokes vector with S', we then have 

n 

S| = X>;^> z = l,. ..,4. (6) 

i=i 

It is useful to derive the expectation value and variance also in the case of S^. Using Eq. (J3]), 



we find 



j=\ fc=i \j=i 



= nAt Si , i 

where we used the fact that V ■ dijrrij k 
Eqs. fl6]) and ©, 



1,...,4, 



(7) 



5jfc. Similarly, for the variance of S^ we find, using 



a\S[)=E(\S[-Stf) = E\ 



-1 2n 



2_^ V^ i 3 



J2 dijd lk E(\Jj - J 3 \[Jt k - J k }) 
j,k=i 



4. 



(8) 



Equations ([3])-([5]) and (J7|) and (jHJ) provide the basic formulas through which we can eval- 
uate the performance of different modulation schemes in the presence of atmospheric seeing. 
They will be related to earlier work in Sect. E] 

When the exposure time At is sufficiently large compared to the coherence time of the 
atmospheric seeing, like in the case of slow modulation cycles, then we can assume that the 
covariance terms E(\J"j — J?j\[^k ~ ^k\) i n Eq. OH]) are negligible for j ^ k. In such case the 
variance of S ; t can be expressed directly as a diagonal quadratic form of the variances of the 
modulated intensity signals, 



<r 2 (sO = E4*V;), 



1,...,4. 



(9) 



i=i 



In Sect. S]we provide a rigorous demonstration of the above statement. However, in many 
cases - certainly when fast cameras with frames rates > 10 Hz are employed - the time 
interval t i+1 — ij (whose inverse conventionally defines the modulation frequency) is of the 
same order or less than the time scale to of the atmospheric seeing. As a result, the inten- 
sity variations induced by seeing in different states of the modulation cycle are in general 
statistically correlated, so one must use the more general expression of Eq. (JHl). 

We note that Eq. (Q also holds in the absence of seeing, in which case a 2 (J^j) evidently 
reduces to just the contribution due to photon noise. Because the signals from different 
camera exposures are always statistically independent, no covariance terms arise in this 
case, and the photon noise only affects the diagonal terms of Eq. (jHJ). This provides a simple 
recipe to include photon noise in the results presented in this paper. Thus, for notational 
convenience, we will simply drop the photon-noise terms in all of the following treatment. 

Finally, we observe that the diagonality of cr 2 (S , J / ), as expressed by Eq. (j2j), is a fundamental 
assumption in the derivation of the optimal form of demodulation matrices presented in [5]. 
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In other words, the usual definition of the efficiency of polarization modulation schemes 
relies on the fact that seeing-induced noise be negligible with respect to photon noise, or 
at least that the seeing-induced covariances of the type appearing in Eq. (JS} be vanishing. 
Since that approach is based on the minimization of the noise of the Stokes measurements 
in the absence of systematic errors, which is expressed by a relation formally identical to 
Eq. Q, there remains an unanswered question: how would the condition of optimality of the 
demodulation matrix change, if the more general expression of Eq. (jSJ) - which takes into 
account the systematic errors due to the seeing - were adopted. This is a subject of research 
in its own right, which we are not going to address further in this paper. 

3. Dual-Beam Polarimetry 

In the case of dual-beam polarimetry, one has two independent sets of n measurements that 
can be combined into 2n new intensity signals of the form ^ i = J?^ ± J^ 6 , where "a" and 
"6" refer to the two beams. For an ideal polarimeter, the modulation vector for the two 
beams satisfy the simple relations m\(t) = m\(t) = 1, and mf(t) = —m^(t) for % = 2,3,4. 
However, in practical cases the two beams are never perfectly balanced. We take into account 
beam imbalance through the detector's gain factor in front of Eq. ([2"lh and introduce new 
modulation vectors for the dual-beam system, m + (t) and m~(t), according to 

Rmf(t) = K a m a (t) ± K b m b (t) , R = K a + K b . (10) 

Through these new vectors, Eqs. ([3])-([5]) can be directly extended to take into account the 
dual-beam redundancy. It is convenient to introduce a new 2n intensity vector with the 
corresponding modulation vector, 



, J~) T , m±(t) = (m + (£), m~(t)f 
We then find, for i — 1, . . . , 2n, 



i 



Jt = E(St) = RAtJ2^ S, , (11) 

o\J±) = E{[S? - Jt\ 2 ) (12) 

4 +At/2 

= k2 J2 ff m t^ + f i) m k (*' + U) E([VSj ■ x(t)][VS k ■ x(t')]) dtdt' . 

j,k=l _ Ai y 2 

Note that we dropped the time shift ti from both arguments of x(t), relying on the station- 
arity of seeing (see the discussion after Eq. [5], and Sect. H]). The index of tj in the argument 
of the modulation functions m • (t) must be modulo n, since it is related to the actual time 
stamp during the modulation cycle. 
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In the usual way, we can associate with M 1 * 1 = (mf-) an "optimal," dual-beam demod- 
ulation matrix, D*. Then, the extensions of Eqs. (JTj) and flEl) to the case of dual-beam 
polarimetry become, respectively, 



2/i 



3 = E4>? 



i = l,. ..,4, 



In 



j,k=i 



4. 



(13) 
(14) 



As before, off-diagonal covariances in Eq. f fl4|) can be neglected only in the limit of exposure 
times much larger than the seeing coherence time (At ^> to). 

In order to illustrate the effect of beam imbalance, we consider the ideal modulation matrix 
for a general modulation scheme with n states, 



/I 



M 



m 12 rai3 m M 



\ 



(15) 



\1 m n2 m n3 m n J 



If we introduce the degree of beam imbalance, p = (n a — Kb)/( K a + K b), the modulation matrix 
M ± appearing in Eq. (TTTj) is accordingly given by 



A 



M 



± 



pm 12 pmi 3 pm u 



(16) 



1 pm n2 pm n3 pm n4 
p rai2 mi 3 mi4 

\p m n2 m n3 m nA J 

where the top and bottom halves correspond respectively to the "+" and "— " linear combi- 
nations of the n measurements from the two beams. 

From the form of M we see immediately that the seeing variations on the intensity and 
those on the polarization parameters of the incoming Stokes vector are decoupled when p = 0. 
Such removal of cross-talk between intensity and polarization is in fact the rationale behind 
dual-beam polarimetry. However, this ideal goal is attained only if the intensity signals of the 
two beams can be balanced (either by camera gain adjustment, or by rescaling of the signals 
during data reduction) with an error which must be better than the target polarimetric 
sensitivity. To clarify this point, let us assume that the gain factors n a and Kb have been 
experimentally determined with some error by the polarization calibration procedure, so 
that K ab = n' ab + Shi' ab , where n' ab is the measured value of K ab . From the definition of the 



dual-beam modulation vector, Eq. ( ITOj) . after rescaling, we have 



<K <K\ / ( n , <*< , <K\ _ 1 /<K <K 



± ^a n , ^b h _ ^a , ^b 

k m,;, = — - m tt ± — m^ , k = — - + — - 
Correspondingly the degree of beam imbalance becomes 

K a K b J I V K a K b J 2 \ «' a K fe 

so that \p\ < (l/2)(|5^|/« / a + \5n b \/K b ) provides a sensible bound on the maximum error 
acceptable in order to achieve a given polarimetric precision in the presence of residual 
intensity-to-polarization cross-talk. 

Since the two beams can in principle always be rescaled post facto, we will assume in 
the following that we are always dealing with perfectly balanced beams. It is understood 
that the possible rescaling of the two beams has to be taken into account for the proper 
determination of the noise on the the combined beams, according to the usual formula 
a 2 (aX + bY) = a 2 a 2 (X) + b 2 a 2 (Y), where a and b are real numbers. This has an effect on 
the derivation of the quantities a 2 (J^^) in Eq. (TT2|) from the noise statistics of the individual 
beams. In the appendix |A] we illustrate this problem for the particular case where the beam 
imbalance is produced by the differential efficiency of a diffraction grating in the p and s 
polarizations. 

4. The behavior of seeing-induced correlations between modulation states 

We must evaluate covariance terms E([^?j — J?j] [J^ — J^]) in Eq. (jHJ) (or Eq. [H]), and study 
their behavior as a function of the modulation frequency. Using the definition (J2J), 

E^-J^-A]) (17) 

4 / +At/2 +At/2 \ 

= k 2 Y^ E I mpit + t^VSp-xit + t^&t I m q (t' + t k )VS q -x(t' + t k )dt' 

p ' q=1 \-At/2 -At/2 ) 

4 +At/2 

= n 2 At 2 J2 d a Sp d S q -^ ft m p (t + 1,) m q {t' + t k )E(x a (t + 1 3 ) x p {t' + t h j) dt dt' , 

p ' 9=1 -At/2 

where in the last line a double summation over the coordinate indexes a, (3 = 1,2 is implicit. 
Because seeing can be considered a stationary random process, the expectation value in 
the last line of Eq. (ITT)) can be written in terms of the two-time correlation matrix, 

T a p{H -t) = E {x a {t) xp(H)) , a, (3 = 1,2, (18) 



and because of the isotropy of the seeing motion, we also hav< 

Tii(t) = T 22 (t) = T(t) , (19a) 

r 12 (t) = r 21 (t), (i9b) 

since the two components of such a motion cannot be distinguishable. In addition, the two 
component motions x\(t) and x 2 (t) are orthogonal, and therefore they can be assumed to be 
independent random processes. Hence, 

E(x a (t)xp(t')) = E(x a (t))E(xp(t')) = , a^P, 

since X\(t) and x 2 (t) are random processes with zero average. Thus the correlation matrix is 
diagonal, and proportional to the unit matrixjj 



T aP {t) = 8 aP Y{t) 



Equation (fT7|) then becomes 



E{[Jj-J 3 \{J k -J h \) (20) 



4 +At/2 



i 2 At 2 J2 VS p • VS q a72 // mp ^ + tj ^ m ^ + tfc ) r ^' ~ t + tk ~ t i) dt dt ' 



p ' g=1 -At/2 

Following the formalism of Sect. [3j for a ra-state modulation scheme in dual-beam configu- 
ration, j and k vary from 1 to In. However, the indexes of tj and tk must be taken modulo 
n, because these refer to the actual steps of the modulation cycle. 

No further simplifications can be made at this point in the case of a continuously mod- 
ulating device. In the remaining part of this section, we will then restrict ourselves to the 
case of stepped modulators, for which the following relation holds 

rriiit + tj) = rriiitj) = m^ , Vt G {-At/2, +At/2) . 

The last equivalence follows directly from the definition of the modulation matrix, Eq. (J3J), 
and so we can simply operate the substitution m p (t + tj) m q {t' + tk) — > mj P mk q in Eq. (120]) . 
Using standard manipulations [1], the double integral in Eq. fTSOj) can then be transformed 



The cross-correlation function of a real stationary random process satisfies the relation T a p (t) = T^ a {—t). 
The additional symmetry constraint provided by Eq. (|19b ) can be viewed as a consequence of the time- 
reversal symmetry of the seeing motion. 

6 Another way to state this result is to consider that since the correlation matrix (|19|) is symmetric, it can 
be diagonalizcd via a similarity transformation involving standard rotation matrices in 0(2). On the other 
hand, because of the isotropy of the seeing motion, there cannot be any preferential direction to attain such 
a diagonal form, and so the correlation matrix for the seeing displacement must always be diagonal. 
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Fig. 1. Two-time correlations of the seeing-displacement normalized ampli- 
tudes, ^ij_fc|(At), plotted against the normalized exposure time, xAt, for var- 
ious values of \j — k\: (variance; thick continuous line); 1 (first neighbor; 
dotted line); 2 (second neighbor; dashed line); 3 (third neighbor; dot-dashed 
line). For this calculation, we used the expression ( 1221) for the auto-correlation 
function T(t), and a duty cycle r = 1 for the camera. 

into a single integral. In fact, noting that tj, — tj = (k — j)At/r, where r is the duty cycle of 
the camera (0 < r < 1), we have 

+At/2 



%,{At) 



At 2 

1 
At 



T(t' - t + sAt/r) dtdt' 



-At/2 
At 



~ J T(r + sAt/r) dr 



(21) 



-At 



For the sake of demonstration, in the following we assume for T(t) a functional dependence 
typical of a Gauss-Markov random process, 

r(t) = r(o) e - x|t| , 



x>o 



22) 



It is known [6] that the Kolmogorov description of atmospheric turbulence leads instead to 
an auto-correlation function of the form e~ x '*' . However, its use in place of Eq. (1221) would 
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Fig. 2. Power spectra of a seeing-like random process. The thin curve repre- 
sents the normalized power spectrum corresponding to the auto-correlation 
function ( 1221) . The thick curve represents a typical modification of this spec- 
trum produced by adaptive optics. 



not change qualitatively the conclusions of this section, since Eq. ( 1221) already contains the 
essential features of the seeing power spectrum that we are going to analyze. Using the 
auto-correlation function ( 122]) . jT s (At) can be integrated analytically. 



#(At) 



r(o) 



r(o) 



>xAt _ !> 



- X At(l+|s|/r) 



x 2 At 2 
2(xAt + e-* A * - 1) 



s^O 



. 



X 2 A£ 2 

We then see that ,%{At) —¥ when At ^> x^ 1 as expected, because of the random nature 
of atmospheric seeing. In particular, for very long exposure times, 3F s (Ab) tends to zero at 
least as (;\A£) -2 for s ^0, whereas ^o(At) ~ (xAt) -1 . For typical atmospheric conditions, 
X^ 1 ~ to ~ 0.01 s, which is of the same order of magnitude of typical exposure times. 
Hence, the terms E([J?j — J^j\[J^k ~ ^k\) must in general be taken into account for a proper 
determination of the measurement errors on the Stokes vector, even for j ^ k. Figure [H 
shows how fast these off-diagonal terms drop as \j — k\ increases. A camera duty cycle with 
r = 1 was assumed for that figure. 
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Fig. 3. Normalized auto-correlation function corresponding to the two power 
spectra of Fig. [2j The thin curve corresponds to Eq. ( |22l) [after normalization 
by T(O)], whereas the thick curve corresponds to the AO-corrected spectrum. 
Note the significant reduction of the time interval within which an efficient 
suppression of the auto-correlation of the seeing is attained in the presence of 
AO correction. 



In order to illustrate the effects of the covariances ( l20|) on the seeing noise on the measured 

V5r V5j, which allows 



Stokes parameters, we introduce the Stokes gradient matrix, G;j 
us to recast Eq. ( 120]) in matrix form, 

Cov jk (At) = E{[Sj - J 3 \\J k - A]) = k 2 A£ 2 (MGM r ) ^_ fc ,(At) . 



(23) 



Here we used Eq. (I2~T]) . and the fact that ^_fc(At) = ,%-j(At), as indicated by the explicit 
functional form of those integrals in practical cases. Consequently, Eq. (jHJ) can also be written 
in matrix form as 

a 2 (S.;) = (DCov(At)D T )^. (24) 

We conclude this section by studying modifications of covariances (120]) arising from the use 
of adaptive optics (AO). The general effect is an important reduction of the low- frequency 
part of the seeing power spectrum. Figure [2] shows an example based on the model of seeing 
correlations described by Eq. ( 1221) . Those curves should be compared qualitatively with the 
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Fig. 4. Same as Fig. [Q but assuming the AO-corrected auto-correlation func- 
tion shown in Fig. |3j Note the significant damping of the covariances for a 
given value of xAt, compared to the case of uncorrected seeing. As for Fig. HJ 
a camera duty cycle r = 1 was assumed for this calculation. 



models of the seeing used in [2]. Here we define the power spectrum S(u) as the Fourier 
transform of the auto-correlation function [I] 



+oo 



S( U ) 



T(t)e- luJt dt = T(0) / 7 (t)e 



Au)t 



dt 



(25) 



where 7 (£) = r(t)/r(0). The thin curve in Fig. [2] represents the power spectrum of the 
uncorrected seeing motion, whereas the thick curve shows a typical modification of this 
spectrum produced by the action of adaptive optics. Note how the high-frequency part of 
the spectrum is not modified by the seeing correction. 

In practical cases, the seeing will be described by an observed power spectrum, S(u), 
rather than by a model auto-correlation function, T(t). From Eq. (|25|) . we then have 



r(t) = — / S(cj)e luJt dio 



(26) 



11 



through which Eq. (12 ip becomes 



+00 



e iu>(k-j)At/r du _ ^ 



- 5( W ; 

7T 



-At 
1 — COS UO At 



'■At 2 



Because S(u>) = S(—u>), it is easy to verify that Eq. (12"7|) implies ^j_k(At) = <%-j(At), as 
we had already derived earlier. 

The application of Eq. (1261) to observed power spectra allows the determination of realistic 
auto-correlation functions of the seeing motion, even in the case of AO-corrected systems. As 
an example, Figure [3] shows the normalized auto-correlation functions corresponding to the 
two power spectra of Fig. [21 Obviously, the thin curve corresponds to the auto-correlation 
function (1221) . Once the true auto-correlation function for an AO-corrected system is known, 
one can use Eqs. (1201) and (12T1) to determine how the modulated intensity covariances are 
modified by the AO correction. Figure |4] shows this effect in the case of the AO-corrected, 
auto-correlation function shown in Fig. [3] (thick curve). Comparing these results to those of 
Fig. HI which correspond to the case of uncorrected seeing motion, we see that the variances 
drop faster in the presence of AO correction. On the other hand, the covariances do not 
converge to zero any faster than in the absence of AO correction. These covariances however 
change sign, so one could in principle adopt modulation schemes and modulation rates such 
that some of the covariances are either vanishing, or even contributing a negative term to 
the expression of a 2 {S[), thus possibly reducing the final error. 

4. A. An illustrative example 

In order to illustrate the effect of seeing on polarization measurements, based on the results 
derived above, we will consider two popular, step-wise modulation schemes: a "Stokes defi- 
nition" scheme and a "balanced" scheme, with ideal modulation matrices given respectively 
by 
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(28) 



Both schemes provide maximum, balanced modulation efficiencies for all Stokes parameters 
(1 for Si, and l/y/3 for S«, with i > 2). The respective modulation matrices in the dual-beam 
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case are derived according to Eq. (TT6|) . We consider here the case of perfect balancing of the 
two beams. 

In the case of the Stokes-definition scheme, by design, only one of the Stokes parameters 
contributes at any time to any given intensity signal combined from the two beams, so that 
all cross-talk terms between different Stokes parameters are eliminated. In addition, the 
seeing- induced Stokes variations enter the Stokes-definition scheme only through terms that 
are diagonal in VSj. In the case of a balanced modulation scheme, instead, one must take 
into account general covariance terms that depend on both ViSj and VS^, for j, k = 2, 3, 4, 
which complicate the expressions of the seeing-induced polarization cross-talk. 

For illustration consider a particular case with VS2 = V1S3 = 0, and VSi parallel to V1S4. 
This example might represent realistic distributions of magnetic fields on the solar surface, 
such as those of magnetic bright points associated with emerging flux, observed near disk 
center. Introducing then the quantities g\ = |VSi| and g± = IVS4I in Eq. ( |24"|) . we find, for 
the Stokes definition scheme, 

a 2 (S[) = \#*Pgl WAt) + f £Li(6-*)#(At)] , (29a) 

a 2 (S' 2 ) = a 2 (S' 3 ) = , (29b) 

a 2 {S' A ) = i n 2 At 2 g 2 \%{At) + #(At)] , (29c) 

while for the balanced modulation scheme, 

o 2 (Si) = J« a At 3 ^[^o(At) + |EL(4-«)^(At)] , (30a) 

a\S' 2 ) = \K 2 At 2 g 2 {%{At) - ST 2 {At) - \{^{At) - ^3 (At)]} , (30b) 

a 2 (S' 3 ) = \ K 2 At 2 g 2 {%(At) - %{At) + |[^(Ai) - ^ 3 (At)}} , (30c) 

a 2 (S' 4 ) = \K 2 At 2 g 2 l%(At) + \Y? s =i^-s)^(At)\ , (30d) 

If we assume that the exposure time, At, is sufficiently large compared to the characteristic 
time of the seeing, t 0) then we can neglect all £^ s (At) terms with s 7^ in the above equations. 
In such case, we see that the two schemes are affected by the same total error on the 
inferred Stokes vector, under identical conditions of camera exposure and time duration of 
the observation. This is because there are 6/4 = 1.5 more modulation cycles for the balanced 
scheme than for the Stokes-definition scheme, during the same time interval. Equations (1291) 
also show that the seeing-induced error in the Stokes-definition scheme only affects the Stokes 
parameters that have non- vanishing gradients at the entrance of the modulator. In particular, 
<?4 only induces an error on 5*4. In contrast, in the balanced scheme, the same error is evenly 
distributed among all of S 2 , S 3 , and £4. 

It is important to remark that the variances expressed by Eqs. (|2"9"1) and (1301 strictly apply 
to Stokes vectors entering the modulator of a polarization-free telescope. In general, instead, 
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typical telescopes' Mueller matrices map the gradient vector (VSi, 0, 0, V1S4) at the entrance 
to the telescope onto a new gradient vector (VS^, VS 2 ) VS3, VS4) at the entrance to the 
modulator, so that all inferred Stokes parameters are affected by seeing-induced errors. In 
such case, there is no evident advantage in adopting the Stokes-definition scheme over a 
balanced modulation scheme. 

5. The effect of seeing during a full observational sequence 

The results of Sect. H] apply to observations that consist of a single modulation cycle, such as 
is commonly employed by imaging polarimeters. The typical mode of operation of slit-based 
spectro-polarimeters instead is to integrate over many modulation cycles for each position of 
the slit in a map. Even at moderate modulation frequencies, the seeing-induced correlations 
in adjacent modulation cycles do not vanish, so we cannot neglect covariance terms between 
different modulation states in different modulation cycles. Hence, we must regard the entire 
elemental observation as a single measurement. In order to do so, we consider the expression 
of the Stokes variances, Eq. (jHJ), and extend it to the case of a series N of modulation cycles. 
Substituting Eq. (120)) . we find 

nN 

j,k=i 

+At/2 



H 



4 nN 

2 



J2 V5 P • VS, Yl di i dik 1 1 m P^ + *i) m $ + tfc ) T W ~ * + tk ~ l o) dt dt ' ■ 

p,q=l j,k=l 

We then introduce the box function 



, 1 , \t\< a/2 
, \t\> a/2 

substitute Eq. (126J1 . and finally operate the changes of variable £, = t + tj and £' — t' + tk- 
We find 

2 4 +2° +2° nN 

° 2 ^ = ^ E vs ? ■ V ^ J duj S ^ J e_i ^ E *** n ^ - *i) m ^) de 



p,q=l 7=1 

— 00 

+00 



/ e 1 ^' Y d * IM£' - h) m q {0 d£' . (31) 
J 1 — 1 



To make this formula applicable to specific modulation schemes and seeing realizations, 



we must evaluate the Fourier transform Hij(u) of the functions 



nN 

H u(0 = 7TT E d * k Ha *^ ~ tfc ) m ^) ■ 



At 
k=i 
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To this purpose, we first note that the elements dik are periodic in k with period n, so we 
can rewrite 

n N-l 



H iM) =mY1 dik J2 Ua ^ ~ tk ~ lnAt / r ) m A0 



At 

k=l 1=0 



'N-l 



x7 Z) dik [ Yl s ^ - tk ~ /nAt / r ) * [ nAi ^) m ^ + tfc )] ) ' ( 32 ) 



At 

fe=l \ l=o 



where in the last equivalence we also used the periodicity of rrijit). 

We take advantage of the convolution theorem of Fourier analysis to derive Hijioo). We note 
that Hij(£) is a real- valued function, so its Fourier transform has the conjugation property 
H*Au) = Hij(-uj). Equation (IHTj) then becomes 



-oo 



K At f 

o-TO = ~^~ E V ^ • VS « J S ^ H ip (u)H* q (u) doo 



2 A +2 



K 2 At 

2n 



p,q=l 

+oo 



J S(u) (H(w)Gflt( W )) ft do; , (33) 



-oo 



where i-T(u;) is the Hermitian conjugate of H(w). Using the symmetry properties of H(u;) 
and G, it is simple to demonstrate that the integrand in Eq. ( 133]) is an even function of oo. 
This allows us to restrict the integration domain to [0,+oo), which is where the observed 
power spectrum is naturally defined. 

The case of a continuously rotating modulator is summarized in the appendix [B], where 
we limit ourselves to providing the Fourier transforms of the functions II a* (^) ^ (£ + £&) that 
appear in Eq. ( 1321) . In the case of a stepped modulator, rrij(£ + tk) — rrikj- For the unit box 
and the "windowed" comb functions in Eq. (1321) . the following Fourier transform pairs hold, 

n a (£-£o) < — ► a e^ sinc(wa/2) , 

Y6(Z-t -la) <-► jye-Mfo+^-Da/gl smc(Nu;a/2) 
f— ' smc(coa/2) 

We thus find, also noting that tk — t\ + (k — l)At/r, 



H.Auo) = Nsmc(ooAt/2) ^A^^Ll e -Mfa+(*-DnAt /a r] y d Mk -i )A t/r (34) 

JK ' v ' ' sinc(nwAt/2r) ^ J K J 

Note that dik is an element of the demodulation matrix corresponding to the extended 
measurement of N cycles. Such matrix contains a factor 1/N with respect to the analogous 
matrix for one cycle. Therefore, we can replace dik with the standard (one-cycle) demodu- 
lation matrix, and drop the factor N in front of Eq. ( 1341) . Because the matrix H(oo) always 
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Fig. 5. Plots of the Stokes errors of Eq. ( 133]) normalized to the incoming in- 
tensity, as a function of the modulation rate. The case shown corresponds 
to the balanced scheme of Eq. ( |28l) in single-beam configuration, with non- 
zero gradients of the Stokes parameters, such that g\ — g± — S\ arcsec -1 and 
92 = 93 = 0.1 pi: continuous curve, cr(S[); dashed curve, cr(S' 2 ); dash-dotted 
curve, o~{S' 3 ); dash-triple- dotted curve, o~{S' 4 ). These plots were calculated for a 
total modulation time T = 10 s, assuming a camera duty cycle r = 1. The thin 
curves are for an observed power spectrum of the seeing, while the thick curves 
show the effects of AO correction (power spectra courtesy of T. Rimmele). 



appears in a product with H'(o;), we can drop all common phase factors from its definition, 
so we can rewrite 

H % M = ^<^/2)^^L± d lk m k3 e -i(*-^/"» , (35) 

where we also used the fact that NnAt/r = T, i.e., the duration of one elemental observation. 
In the case of dual-beam polarimetry, the summation in Eq. (|35|) is extended to 2n, while at 
the same time (k — 1) in the exponential must be taken modulo n. 

We note that H ip (u)H* q (u) — > S ip 5i q for vanishing At. Therefore, only Gu = IVSjI 2 
contributes to o~ 2 (Sl) in Eq. (I3"3"|) for At — > 0. This is in complete agreement with the result 
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Fig. 6. Same as Fig. [5j but for the dual-beam configuration with perfectly 
balanced beams. 



that can be derived from Eqs. ( )23|) and (|24|) under the same limit conditions. The results 
presented in [HE] correspond evidently to the diagonal case p = q. 

Figure [5] shows an example of seeing-induced errors on the measurement of the Stokes 
vector for a balanced modulation scheme in single-beam configuration, for a total modulation 
time of 10 s, and with gradients g\ — g^ — Si arcsec -1 and gi = g$ = 0.1 g\. Measured 
power spectra of seeing-induced image motions both with and without AO correction were 
used to produce these plots (T. Rimmele, private communication). This figure shows the 
clear benefit of AO correction in terms of a reduction of the seeing-induced errors by more 
than an order of magnitude for S\ and S4. For modulation periods larger than the seeing 
correlation time, i.e., for long exposures, the variances of S2, S3, and S4 are dominated by 
cross-talk from gradients in Si. When the modulation frequency is increased to the point 
that seeing-induced displacements are practically frozen for the duration of a modulation 
cycle, the cross-talk terms in each error become negligible compared to the diagonal terms. 
The use of AO implies shorter correlation times of the seeing-induced displacements, and 
hence a higher modulation frequency is required to satisfy this condition. 

For the dual-beam case (Fig. |6]) with perfectly balanced beams, there is no cross-talk from 
Si to S 2 , S 3 , and S 4 (see Sect. |3]). This fact is illustrated by the drop in the error curves for 
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Fig. 7. Same as Fig. [6j but for the case of a stepped modulator (with 8 states) 
consisting of a waveplate with 150° retardance, identical to the one considered 
in pp. The same demodulation scheme of [I] was adopted for this calculation. 
These plots are for a dual-beam configuration with perfectly balanced beams. 
The gradients of the Stokes parameters are g^ = gs = S± arcsec -1 and g± = 
#4 = 0. The red curves (see on-line version of the figure) show the predicted 
cross-talk level neglecting terms of the form #2<73) which are missing in the 
treatment of [HE]. 



5*2, 5*3, and S4 with respect to the single-beam case of Fig. [51 However, the cross-talk among 
5 2 , 5*3, and S4 remains. Analogously to the single-beam case, at high modulation frequencies, 
the seeing-induced error on each element of the Stokes vector tends to the contribution from 
the diagonal term only. Comparing the plots of Figures [5] and El we see that for very high 
modulation frequencies (> 100 Hz without AO, or > 1 kHz with AO) the performance of 
single- and dual-beam modulation schemes become comparable, and the only benefit of a 
dual-beam polarimeter in such case is the redundancy of polarimetric information, leading 
to a reduction of the photon noise by a factor of a/2 with respect to the single-beam case. 

Our formalism produces results that agree qualitatively with those of [TJE] , although with 
some notable quantitative differences. Figure [7] shows the polarization cross-talk in a dual- 
beam configuration, calculated for a stepped modulator consisting of a rotating waveplate 
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with 150° retardance, identical to the one considered in the study by [T]. The case considered 
in the figure corresponds to the presence of spatial gradients in S2 and S3 only. The red curves 
show the cross-talk that is derived by neglecting terms depending on VS2 • VS3, which are 
missing in the treatment of [HE] (see comment at the end of the next section). 

6. Discussion 

In this paper, we have approached the determination of seeing- induced cross-talk noise from 
a statistical point of view. Equations (IS])-© and ([7|) and (jSJ) form the basis of our derivation 
(see Sect. EJ), so it is important to compare those results with previous work [HE]. 

In [HE] , the "variances" there defined depend explicitly on the total observation time and 
the modulation frequency, through the integral of the product of the seeing power spectrum 
with a real function of frequency that depends implicitly on both (cf. Eq. [15] of [TJ, and 
points 1-5 of Sect. 2 of [2j). In our formalism, the total integration time does not appear 
explicitly because Eq. (jSJ) represents the variance on a single measurement of J"i. However, 
if a number N of such measurements are made, and those measurements can be considered 
statistically uncorrelated, then the variance on the average signal J^ is reduced by a factor 
N. In such case, for a fixed modulation frequency, but increasing the total observation time, 
we expect the same qualitative scaling of variance with the integration time as derived in 
previous work. Secondly, the dependence on the modulation frequency and the seeing power 
spectrum is also apparent within our formalism - already in the case of a single measurement 
- when we consider Eq. (J5J). 

Incidentally, we note how the power spectrum of the seeing appears most naturally in 
our approach as the Fourier transform of the two-time correlation function of the seeing 
displacement vector (see Sect. 0j), whereas in previous work it is identified instead with the 
modulus square of the Fourier transform of the seeing displacement J3 The correspondence 
between these two approaches to the definition of the power spectrum of a stationary random 
process is clearly described in pEJ. 

The importance of the covariances E([J^j — J'^^y. — ^]) in typical cases has fundamen- 
tal implications for the concept of polarimetric measurements. Through our analysis we are 
able to quantify the significance of seeing-induced correlations between measurements cor- 
responding to different modulation steps, and how these correlations decay for decreasing 
modulation frequencies. Based on those results, it must be expected that seeing-induced 
correlations between different modulation states typically extend beyond the time interval 



We observe that the ordinary Fourier transform is not well defined for the seeing displacement vector, 
because the associated random process is not limited in time. One can get around this problem by defining 
a generalized Fourier transform of the seeing, as the limit of ordinary Fourier transforms of finite samples of 
the seeing for increasing duration of those samples [I] . 
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of one modulation cycle. In other words, the Stokes vector measurements corresponding to 
different modulation cycles during an elemental observation are in general statistically cor- 
related. Under this condition, we must expect that the variance of the average signal J?i will 
obey the ordinary scaling law by the total number N of cycles of the elemental observa- 
tion only approximately. Formally, one should consider instead such elemental observation 
as a single measurement that is realized through the totality of the N modulation cycles, 
and thus characterized by a corresponding nN x 4 modulation matrix. The variance of this 
measurement can then be determined through the usual equations, adopting such extended 
definition of the modulation and demodulation matrices (see Sect. [5]). A similar approach is 
taken also in the earlier work [HE], where these correlations are made manifest by performing 
a spectral analysis of the demodulated signal over the entire time interval of the elemental 
observation. 

Our formalism reveals however that in [HE] the covariance terms corresponding to con- 
tributions proportional to VSj ■ VSj, with i ^ j (see Eq. [20]), are not accounted for. In 
fact, in [T] the variances are written down as Eq. (12) of that paper directly from Eq. (11), 
which are explicitly of diagonal form. In other words, looking at our Eqs. ( 1231) and ( 1241) . 
previous work has computed the seeing- induced noise always under the assumption that the 
gradient matrix G was diagonal. The additional off-diagonal components would arise in the 
development of [HE] when properly evaluating the expectation value of (O r — O r ) 2 in the 
notation of [TJ. Equations ( 1231) and ( 1241 also clarify the vanishing of the cross-talk terms in 
the limit of large modulation rates (x^t — > 0). Because all the integrals ^L_^|(At) tend to 
the same value in such limit, the vector of the Stokes variances becomes simply proportional 
to the diagonal of the gradient matrix, G. 
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A. Beam imbalance introduced by a polarizing grating 

We consider a beam of unpolarized light incident on a diffraction grating. Typically the 
diffraction efficiency of a grating is different in the p and s directions, so the diffracted 
beam consists of a mixture of two orthogonally polarized beams with generally different 
intensities S p and S 8 . Let us consider the case in which these two beams get separated by 
a perfectly polarizing beam-splitter placed after the grating. The emerging beams will then 
have intensity S p and S a , respectively. 

Let us consider the combined signal of intensity S so defined, 
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where k p and k s are the two gain factors applied to the detector to balance the two beams. 
The ratio k p /k s then corresponds to the ratio of the two orthogonal efficiencies of the grating, 
e p and e s , according to 

Kp/^s = e s /e P • (36) 

The photon noise associated with the combined signal S is given by 

- (t) v < s >> + (t) v < s -> • 

Correspondingly, the relative error is given by 

a 2 (S) fn p S p \ 2 a 2 (S p ) , f k s S s \ 2 a 2 (S s ) 



C2 \ 77. q / C2 \ 77. q / C2 



Evidently, for balanced beams 
and therefore 



ftp'Jp — K s>->s — 2 ' (37) 



a 2 (S) 1 



a 2 (S p ) a 2 (5 s 



C2 ' C2 



(38) 



S 2 4 

If we express the signals in terms of photon flux, assuming Poisson's statistics for the photon 
counts, and indicating with r N the read-out noise of the camera, we can rewrite Eq. (158]) as 

a 2 (S) 1 fN p + r 2 N N s + r 2 N * 



S 2 4 V N 2 N 2 



If we indicate with N in the photon count before the grating, evidently N p = ^e p N in and 



N s = |e s iV"i n , so we find 



v 2 (S) 1 e p + e r 2 N e 2 + e s 2 

S 2 2N- m e p e s + N 2 n e 2 p e 2 

1 F r 2 £ 2 , 2 

= — — + ^- p 8 (40) 

iV- e e JV- 2 e 2 e 2 ' l ; 

Jv in c p c s Jv ln c p c s 

where we also introduced the grating's average efficiency e — |(e p + e s ). 

We note that the applicability of the above treatment in practical cases relies on the 
previous knowledge of the grating efficiencies, e p and e s . These must be determined as a 
function of wavelength from flat-fielding, ideally using unpolarized radiation. If the incident 
beam is instead weakly polarized (as it is typically the case for sunlight flat-fields), the 
above treatment is still valid, but Eq. ( 136]) can only provide an approximate estimate of the 
appropriate gain factors, k s and /%,, if the measured efficiencies of the grating are adopted 
in that equation. 
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B. Rotating modulator 

We consider the case of a retarding device, which is continuously rotating with angular 
frequency Q, and which at 0° position is described by the Mueller matrix 

/l \ 

fJ,22 /^23 ^24 

fJ-32 /^33 ^34 

\0 /i42 /i43 H44/ 

We note that a full modulation cycle corresponds to only half rotation of the modulator, 
because of the characteristic 180°-periodicity of polarization modulation. 

In order to determine the cross-talk terms f|33|) for such a device, we need to compute the 
Fourier transforms of the functions (see Eq. |32j ) 

z jk(0 = IWO mj(£ + t k ) , j = 1, . . . , 4 , fc = 1, . . . , n , 

where n is the number of modulation states (camera exposures) in the modulation cycle. 
These Fourier transforms are given by the following expressions, 

Z lk (ou) = At sinc(u;At/2) , 

Z 2 k{u) = - (/i22 + ^33) At sinc(wAt/2) 

+\ [/i 22 - ^ + i(M23 + /i 32 )] e 41 ^^-^" 1 )/"] At sinc([u; - 4Q]At/2) 

+\ [m ~ ^3 - i(A*23 + M] e" 41 ^^-^- 1 )/"] At sinc([u; + 4Q]At/2) , 

^3fc(w) = - (/i 2 3 - A*32) At sinc(u;At/2) 

-\ W - W + K/^23 + /i 32 )] e 41 ^ 1 ^^- 1 )/"! At sinc([^ - 40]At/2) 
+ \ [/i 22 - ^33 - i(M23 + /i 32 )] e" 41 ^^-^- 1 )/"! At sinc([u; + 40]At/2) , 

Z ik {uj) = - (/i 24 + i/i 3 4)e 2i[ ™ 1+7r(fc - 1)/n] At sinc([w - 2Q]At/2) 

+ - (/i 24 - i/i 34 )e' 2i[ntl+ " (fc - 1)/nl At sinc([w + 2Q]At/2) , 

where At is the exposure time for each modulation state, and ti is the time at which the 
modulator is found in the 0° position. 
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